Solving Einstein's Equations With Dual Coordinate Frames 
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A method is introduced for solving Einstein's equations using two distinct coordinate systems. 
The coordinate basis vectors associated with one system are used to project out components of the 
metric and other fields, in analogy with the way fields are projected onto an orthonormal tetrad 
basis. These field components are then determined as functions of a second independent coordinate 
system. The transformation to the second coordinate system can be thought of as a mapping 
from the original "inertial" coordinate system to the computational domain. This dual-coordinate 
method is used to perform stable numerical evolutions of a black-hole spacetime using the generalized 
harmonic form of Einstein's equations in coordinates that rotate with respect to the inertial frame 
at infinity; such evolutions are found to be generically unstable using a single rotating coordinate 
frame. The dual-coordinate method is also used here to evolve binary black-hole spacetimes for 
several orbits. The great flexibility of this method allows comoving coordinates to be adjusted with 
a feedback control system that keeps the excision boundaries of the holes within their respective 
apparent horizons. 



I. INTRODUCTION 

This paper introduces a new method of solving Ein- 
stein's equations using two distinct coordinate systems. 
Tensors like the metric are represented by their compo- 
nents in the coordinate basis of the first coordinate sys- 
tem. Einstein's equations are then used to determine 
these tensor components as functions of the second co- 
ordinate system. The mapping between the first and 
second coordinate systems is quite arbitrary and can be 
chosen dynamically. This freedom allows us to adapt the 
second coordinate system continuously, for example to 
track the motion of the individual black holes in a bi- 
nary. The first coordinate system plays much the same 
role as the orthonormal tetrad used in some formulations 
of the Einstein equations [3, 0, El • The transformation 
to the second coordinate system can be thought of as a 
mapping between the original "inertial" coordinates and 
the computational domain, and plays much the same role 
as the "grid velocity" sometimes used in numerical hy- 
drodynamics [U. The bulk of this paper gives a careful 
description of this new method as it applies to the gener- 
alized harmonic form of the Einstein equations [j| , plus 
a set of numerical tests that demonstrate its usefulness. 

The flexibility to choose a coordinate frame adapted 
to a particular physical problem is often used to simplify 
solutions of the Einstein equations. For example, solu- 
tions having some symmetry {e.g., time- independence) 
are much simpler when expressed in coordinates that re- 
spect that symmetry. It has long been expected that bi- 
nary black-hole spacetimes would most naturally be rep- 
resented in a coordinate frame that corotates with the or- 
bit of the holes; this frame should be advantageous for nu- 
merical simulations since it would make the fields nearly 
time-independent during the inspiral. Unfortunately our 
attempts to model binary black-hole spacetimes in the 
conventional single-coordinate framework have all failed. 



These failures reveal several serious and interesting prob- 
lems with the single-frame approach, and motivated us 
to develop the new methods described in this paper. Be- 
fore we turn to our discussion of the new dual-frame ap- 
proach, however, we believe it is useful to describe those 
problems with the standard single-frame methods. 

The groundbreaking binary black-hole evolutions of 
Pretorius @ use a single coordinate frame that asymp- 
totically approaches the inertial frame at infinity. The 
singular interior regions of the holes are excised from 
the computational domain. In this approach the black 
holes move across the coordinate grid and the solution 
is time dependent on the orbital timescale. We have at- 
tempted to implement a similar scheme using our spec- 
tral evolution code, but encountered two significant prob- 
lems: First, the black-hole "excision" boundaries must 
be changed from time to time to track the motions of 
the black-hole horizons. This requires finding a way of 
adding (and deleting) grid points to (and from) the com- 
putational domain in a way that keeps the distance be- 
tween the horizon and the excision boundary more or less 
fixed, and then — the difficult part — to find sufficiently 
smooth and stable ways to fill {e.g., by extrapolation) 
the new grid points with appropriate values for all the 
dynamical fields. This problem has been solved suc- 
cessfully for lower-order finite-difference numerical meth- 
ods SSBBO- The difficulty for the spectral-method 
implementation seems to be the inability to produce suf- 
ficiently smooth and well-behaved extrapolations of the 
needed dynamical field quantities to points lying beyond 
the current computational domain. Our attempts to de- 
velop such techniques have not been successful. 

There is also a second fundamental problem with 
the conventional moving-excision method: although the 
black-hole horizons move through the computational 
grid, the excision boundaries are fixed except when they 
are moved to update the grid. To understand why this 



can cause trouble, consider first an excision surface lo- 
cated just inside the horizon, but moving continuously 
along with the horizon. In this case, the excision bound- 
ary is a spacelike surface that lies in the future of the 
computational domain, so (for a causal hyperbolic repre- 
sentation of the Einstein equations) boundary conditions 
are not needed there, and the spacetime region inside 
can be excised. Now consider the conventional moving- 
excision method, where the boundaries are fixed in the 
coordinate grid during a time step, and therefore move 
relative to the horizons. At the trailing edge of the hori- 
zon, the excision surface is moving out of the black hole 
superluminally. Locally this surface is spacelike and lies 
in the future of the computational domain, so boundary 
conditions are not needed there. But at the leading edge, 
the excision surface is falling deeper into the black hole; if 
it falls quickly enough, some part of this surface may be- 
come timelike, so boundary conditions will be required at 
these points. If appropriate boundary conditions are not 
supplied, the evolution problem becomes ill-posed and 
the black-hole excision paradigm breaks down. We moni- 
tor the characteristic speeds at excision boundaries in our 
code, and have verified that this breakdown does occur 
in our simulations when a black hole moves too quickly 
through the coordinate grid. This problem can be ame- 
liorated by moving the excision boundaries deeper into 
the black- hole interiors (where with an appropriate choice 
of gauge, light cones are "tipped" further toward the sin- 
gularity); we presume this is how the successful finite- 
difference codes control this problem. But we were not 
successful in curing this problem in our spectral evolution 
code. Another successful approach to moving black holes 
through the computational grid is the moving-puncture 
method, cf. Refs. [U El ES El OS Ell - which can be 
thought of as shrinking the internal excision boundaries 
to isolated points that lie between the grid points (and so 
are simply ignored). Unfortunately the dynamical fields 
are not smooth at the puncture points, so the moving 
puncture approach may not be extendable to spectral 
methods. 

Since we were unsuccessful in moving black holes 
through our spectral computational domain using the 
conventional method, we were led to pursue an approach 
using comoving coordinates where the black holes re- 
main at rest. Unfortunately, the generalized harmonic 
evolution system [f| exhibits severe instabilities (in the 
sense that the constraints grow without bound) when 
initial data are evolved using coordinates that rotate 
with respect to the inertial frame at infinity. 1 These in- 



1 This instability appears to be a feature of the generalized har- 
monic system. No instability was seen in earlier evolutions of 
black holes in rotating frames using other forms of the Einstein 
system [17|. and we do not find such an instability using our 
code with the KST evolution system [lSj . In the GH system 
the full 4-metric, including the lapse and shift, are evolved as 
dynamical fields, while in the KST system only the 3-metric is 
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FIG. 1: Constraint violations in evolutions of a Schwarzschild 
black hole using the generalized harmonic system and coordi- 
nates that rotate with angular velocity fi. 

stabilities occur even for simple time-independent cases 
like Minkowski or Schwarzschild spacetime. For exam- 
ple, Fig. [1] shows the constraint violations for evolutions 
of Schwarzschild in a frame rotating with angular ve- 
locity fl. These evolutions were performed on a com- 
putational domain that is a little larger than the size 
needed for binary black-hole evolutions, with outer ra- 
dius i?max = 1000M; see Sec. HH for details. Binary 
black-hole evolutions must be performed stably and ac- 
curately for times of order 1000M to model all of the 
interesting inspiral and merger dynamics. However, the 
rotating frame evolution shown in Fig. [1] with £1 = 0.2/M 
(comparable to the maximum angular velocity achieved 
by a black-hole binary just before merger) is unstable in 
the unacceptably short time of about 10 _4 M. We do 
not fully understand the instability in evolving asymp- 
totically flat spacetimes in rotating coordinates using the 
generalized harmonic evolution system. We have reduced 
the severity of this problem with methods discussed in the 
Appendix, but we were not able to solve completely the 
problem of evolving with a single corotating coordinate 
frame. 

Another fundamental limitation on using rotating co- 
ordinates in asymptotically-flat spacetimes is a numeri- 
cal resolution issue: some components of the four- metric 
grow asymptotically like p 2 £l 2 for large values of the 
cylindrical radial coordinate p (to leading order), while 
other components approach M/r, for large values of the 
spherical radial coordinate r. On computational do- 



evolved. Since only the shift grows near spatial infinity in a ro- 
tating frame, the non-dynamical treatment of the shift in the 
KST system may account for its rotating frame stability. 
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mains that extend to r w 1000M and angular veloci- 
ties « 0.2/M (appropriate for binary black-hole space- 
times), these values differ by almost eight orders of mag- 
nitude, so that the dynamical range needed to resolve the 
various field components with sufficient accuracy is diffi- 
cult to achieve using double precision numerical methods. 
This difficulty goes away if one uses the inertial frame 
components of the various fields [6j . 

The dual-frame approach proposed here corrects or cir- 
cumvents all of the problems associated with the more 
traditional approaches that we know about. The re- 
mainder of this paper is organized as follows. The dual- 
coordinate frame method is presented more completely 
in Sec. [ill This new method is applied in Sec. 11111 to the 
case where the two coordinate frames rotate uniformly 
with respect to each other. We test this new rotating- 
coordinate method by evolving Schwarzschild initial data 
in coordinates that rotate with respect to infinity, find- 
ing that the evolutions of Fig. Q] become stable and con- 
vergent. A number of technical numerical issues associ- 
ated with the dual-coordinate method (e.g., how to con- 
struct the appropriate spectral filters) are also discussed 
in Sec. IIII1 Section [TVl further tests this new method by 
evolving a binary black-hole spacetime in uniformly ro- 
tating coordinates. This test fails when the fixed angular 
velocity of the rotating frame no longer tracks the mo- 
tion of the black holes to sufficient accuracy. In Sec. [V] 
we present a more sophisticated application of the dual- 
coordinate method by constructing a flexible coordinate 
map whose parameters are adjusted through a feedback 
control system that keeps the black-hole horizons cen- 
tered on the exesion surfaces. Using these new horizon- 
tracking coordinates, we are able to evolve a binary black- 
hole spacetime in a stable and convergent manner for 
about 4.6 orbits. Our results are summarized, and a dis- 
cussion of possible directions for further development are 
given in Sec. IVIl Finally, we review in the Appendix 
our attempts to make rotating single-coordinate frame 
evolutions stable for the generalized harmonic evolution 
system. 



II. DUAL COORDINATE FRAMES 

Consider a first-order representation of the Einstein 
evolution system, such as our formulation of the general- 
ized harmonic system Q . The evolution equations for the 
dynamical fields u a for such systems can be represented 
abstractly as 



(1) 



where A ka p and F a may depend on u a but not its 
derivatives. We use Greek indices a, j3, ... to la- 
bel the various dynamical fields, and Latin indices z, 
j, k, ... to label the spatial coordinates. The bars on 
the various indices in Eq. fT} indicate that this sys- 
tem evolves the coordinate components of the collec- 
tion of dynamical fields, w Q , as functions of coordinates 



x a = (t,x l ). In the generalized harmonic system, for 
example, u & = {tfj ah , LI a6 , % a5 } where ^- al is the four- 
metric, U ab = -t c d 5 ip a - b , = d- k t(j a - b , and t c is the 
unit timelike normal to the t ^constant hypersurfaces. 
We use Latin indices from the beginning of the alphabet, 
a, 6, c, ... to label spacetime coordinates. 

Einstein's equations are covariant, so it is straightfor- 
ward to transform any representation of those equations 
from one coordinate frame to another. The standard 
method of transformation changes both the coordinates, 
x a — » x a , and the components of the dynamical fields, 
u a — > u a , using the appropriate transformation rules. 
For the generalized harmonic system, for example, the 
field components transform according to the rules: 
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Here we propose a new way of solving the Einstein 
equations that solves both the rotating frame instabil- 
ity problem and the moving-excision-boundary problems 
discussed in Sec. HI We introduce a second coordinate 
system x a that is (in principle) completely independent 
of the first x a . We think of the first coordinate system 
x a as inertial coordinates that do not rotate with respect 
to the asymptotic inertial frame at spatial infinity. We 
use these inertial coordinate bases, <9 a and dx a , to con- 
struct the components of the various dynamical fields: 
u a . These inertial-frame components are well behaved 
near spatial infinity, and the numerical dynamical range 
needed to represent them is significantly reduced. We 
think of the second set of coordinates, x a , as comoving 
coordinates chosen to minimize the time dependence of 
the dynamical fields in some way. We solve the evolution 
system for the inertial-frame dynamical-field components 
u a as functions of the comoving coordinates x a . For sim- 
plicity we consider only the case where the two coordi- 
nate systems have the same time slicing: t — i. In this 
case the system of evolution equations for u a in terms 
of x a is just Eq. JT]) with the straightforward change of 
independent variables: 



d t u & 



dx^ A 



divP = F a . 



(5) 



Here dx % j&t and dx l /dx k are to be determined as func- 
tions of t and x l from the transformation that relates the 
two coordinate systems: x l = x l (t,x k ). 

The dual-coordinate evolution system, Eq. ([5]), has 
many properties in common with the original, Eq. ([1]) . Its 
solutions are the same (assuming the coordinate transfor- 
mation is sufficiently smooth), just expressed in terms of 
the new coordinates x a : u a (x a ) — u a [x b (x a )\. In addi- 
tion the characteristic fields are exactly the same for the 
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two systems, and the characteristic speeds differ only by 
v — v + nidx 1 1 1 dt, where v represents any of the char- 
acteristic speeds of Eq. ([T]) and v is the corresponding 
characteristic speed of Eq. ([5]) . Here rii is the appropriate 
normal one-form used to define the characteristic speeds 
(e.g., the outward-directed normal to the boundary sur- 
face) . Thus boundary conditions for the two systems are 
transformed in the obvious way. Since the characteristic 
speeds may be different, however, it is possible that the 
list of characteristic fields needing boundary conditions 
at a particular point may change. 

III. ROTATING COORDINATE FRAMES 

In this section we present a simple test of the dual- 
coordinate-frame idea described in Sec. |TTJ We evolve 
asymptotically flat spacetimes, like the Minkowski and 
Schwarzschild geometries, using two frames: an inertial 
coordinate frame that is asymptotically Cartesian at in- 
finity, and a comoving coordinate frame that rotates uni- 
formly with respect to the inertial frame. 

Before we discuss those tests, however, it is appropriate 
to describe the numerical methods that we use. These 
tests are done with the generalized harmonic evolution 
system as described in Ref. Q using spectral numerical 
methods as described, for example, in Refs. [HI, [2(|. The 
components of the various fields in these tests are ex- 
panded in terms of scalar spherical harmonics of the an- 
gular coordinates (with I < 11) and Chebyshev polyno- 
mials of the radial coordinate log r through order N r — 1 
(with N r = 15 for the tests in Fig. []}. The computa- 
tional domain used in these tests (as well as those shown 
in Fig. Q} consists of a set of eight nested spherical shells 
with boundaries located at the radii 1.8, 8, 35, 70, 140, 
229, 374, 612, and lOOOAf respectively. We measure con- 
straint violations in these tests with a quantity ||C|| de- 
fined as the ratio of the L 2 norm of all the constraint 
fields of the generalized harmonic system, divided by the 
L 2 norm of the spatial gradients of the various dynam- 
ical fields dkU a (see Eq. [71] of Ref. [5j). This quantity 
||C|| vanishes whenever the constraints are satisfied, and 
the normalization is chosen so that ||C|| becomes of order 
unity when constraint violations dominate the solution. 

In 3+1 formulations of general relativity, the gauge 
freedom in the theory is usually parameterized by a lapse 
function and a shift vector that are freely specifiable. In 
contrast, the gauge freedom in generalized harmonic for- 
mulations of Einstein's equations is represented by four 
freely-specifiable gauge source functions H a (see Ref. @) 
which determine the evolution of the lapse and shift. 
Once H a has been chosen, the harmonic constraint equa- 
tion 

= c a = r a + H a (6) 

must be satisfied. Here r a ee tp bc r a i, c is a trace of the 
Christoffel symbols. For single-coordinate-frame numer- 
ical tests that aim to reproduce a known analytic time- 



independent solution, we typically choose H a so that the 
constraint, Eq. ([B]), is satisfied initially, and we fix this H a 
for all time; this was done for example in the evolution 
shown in Fig. [T] 

For dual-frame evolutions we must be more careful 
with the gauge source functions, particularly because H a 
does not transform like a tensor. Because we have had 
considerable success (in nonrotating frames) choosing H a 
to be time-independent, we take a similar approach here. 
We first define a new quantity H a that has the following 
two properties: 1) H a transforms like a tensor, and 2) 
in inertial coordinates H„ = H a . As in the single-frame 
case, we choose H a so that the constraint Eq. fl6|) is sat- 
isfied initially, but now we demand that H a is constant 
in the moving frame, i.e., that dtH a = 0. 

Our first test of the dual-coordinate-frame idea con- 
sists of evolving Schwarzschild initial data with uniformly 
rotating coordinates. In particular we repeat the most 
unstable evolution shown in Fig. [1] as a dual-coordinatc- 
frame evolution. The inertial coordinates (t, x, y, z) for 
this test are the standard asymptotically Cartesian coor- 
dinates associated with the Kerr-Schild representation of 
the Schwarzschild geometry, while the "comoving" coor- 
dinates (t, x, y, z) rotate uniformly with respect to these 
inertial coordinates: 

t = i, (7) 

x = ircos(nt) + y sin(fit), (8) 

y = — xsin(fit) + ycos(f2i ), (9) 

z = z. (10) 

(The angular velocity, fi, of these co- moving coordi- 
nates can be chosen arbitrarily.) We solve the dual- 
coordinate form of the evolution equations, Eq. (JSJ), for 
these Schwarzschild initial data with £1 = 0.2/M. This 
value of f2 is chosen because it corresponds roughly to the 
orbital angular velocity of a binary black-hole system at 
the time of merger. Figure shows the constraint viola- 
tions for these evolutions for several values of the radial 
resolution parameter N r . These dual-coordinate- frame 
evolutions persist for about 20M, five orders of magni- 
tude longer than the analogous single coordinate evolu- 
tion shown in Fig. [1] Unfortunately these evolutions are 
still unstable on a timescale that is several orders of mag- 
nitude shorter than needed, and this instability shows 
some signs of non-convergence. 

The dual-coordinate-frame representation of tensor 
fields has some unpleasant features that are ultimately 
responsible for the instability seen in Fig. [5] If the 
transformation between the two coordinate frames x l = 
x l (i,x^) is time-dependent (as is the case for uniform 
rotation), then clearly the transformation relating the 
tensor bases, e.g., di — dx J /dx T dj, will also be time- 
dependent. So even if the comoving components of a 
tensor ipab{x J ) are time-independent, as is the case for 
the Schwarzschild example above, then the inertial frame 
components T/> g g(f, ir- 7 ) that we are evolving will change 
with time. Consider for example a purely radial vector 
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FIG. 2: Constraint violations in dual-coordinate- frame evolu- 
tions of Schwarzschild, with comoving coordinates that rotate 
uniformly with angular velocity Q = 0.2/M. Angular filtering 
(see text) is performed on the inertial frame components. 
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FIG. 3: Time-dependence of the inertial frame components 
of two cross sections (equatorial and polar) of a vector field 
expressed as functions of rotating coordinates. At t = the 
vector field is radial, but the d x and dy basis vectors rotate 
with angular velocity Q. The vectors on the right side of the 
polar cross section figure at t — n/(2Q) are pointed into the 
plane of the figure, while those on the left point out. 



field v l (x J ) that is time-independent when represented in 
a single coordinate frame that rotates with angular ve- 
locity f2. Figure [3] illustrates the time-dependence of the 
inertial frame components v l (x 3 ) of this vector field. Two 
cross sections of this field (one equatorial and one polar) 
are shown at time t = 0, when the rotation map is just 
the identity, and at times t = 7r/(2f2) and t = 7r/f2, when 
each inertial frame basis vector <9j is rotated by tt/2 and 
7T about the z axis respectively. 

In addition to making tensor field components time- 
dependent, the dual-coordinate-frame representation can 
also mix the tensor-spherical-harmonic components of 




27I/Q 3n/Q 4n/Q 
t 

FIG. 4: Tensor-spherical-harmonic components of the 
Schwarzschild spatial three-metric expressed in rotating co- 
ordinates. Dashed curves show the components expressed in 
rotating-frame tensor spherical harmonics (only the £ — 
component is non-zero); solid curves show the inertial- frame 
tensor-spherical-harmonic components. 



tensors. In Fig. [3] for example, the radial vector field 
is pure t = Q when represented in rotating-frame vector 
spherical harmonics, but it becomes a time-dependent 
mixture of I = 0, 1, and 2 when represented in inertial- 
frame vector spherical harmonics. Figure [3] illustrates 
the time dependence of the tensor-spherical-harmonic 
components of the Schwarzschild spatial three-metric 
expressed in rotating coordinates. The dashed curves 
are the rotating-frame tensor-spherical-harmonic compo- 
nents, while the solid curves are the inertial-frame tensor- 
spherical-harmonic components. In the rotating frame 
the Schwarzschild spatial three- metric is pure i = 0, 
while the inertial frame components are a time-dependent 
mixture of £ = 0, 1, 2, 3, and 4. The periodicity of the 
tcnsor-spherical-harmonic components in Fig. [4] is deter- 
mined by the rotation period of the coordinates 2ir/Cl. 

This mixing of the inertial-frame tensor-spherical- 
harmonic components is the cause of the instability seen 
in Fig. [2] When using spectral methods in a spherical 
computational domain, we expand the Cartesian compo- 
nents of tensors in scalar spherical harmonics, and we 
evolve these Cartesian components. In this case we find 
it necessary periodically to set to zero the highest-order 
iensor-spherical-harmonic coefficients of all tensors in or- 
der to avoid numerical instabilities [TH, [2l[ . We call this 
operation "filtering" , since it is similar to the filtering 
operations commonly used in spectral methods to avoid 
instabilities caused by aliasing. Filtering is necessary be- 
cause spatial differentiation couples different values of 
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the I and m indices in the scaZar-spherical-harmonic ex- 
pansions of the Cartesian components of tensors. Con- 
sequently any truncated series expansion will result in 
incorrect evolution equations for the highest angular 
modes, and these errors often lead to non-convergent 
instabilities. However, it is possible to truncate the 
iensor-spherical-harmonic expansion of a tensor at a fi- 
nite value of t in a self-consistent way, because the spa- 
tial gradient of a tensor spherical harmonic is also a ten- 
sor spherical harmonic with the same spherical harmonic 
index, but one higher tensor rank. Thus we perform 
filtering by transforming each tensor from a Cartesian- 
component basis to a tensor-spherical-harmonic basis, ze- 
roing the tensor-spherical-harmonic coefficients for values 
of £ larger than those kept in our expansion, and trans- 
forming back to our Cartesian-component basis. This 
filter cures the angular instability problems associated 
with the evolutions of single-coordinate frame spherical- 
harmonic representations of tensors 0, [2l[ . 

This filtering algorithm must be modified for dual- 
coordinate-frame methods because there is now a more 
complicated relationship between the coordinates and the 
basis vectors used to represent tensors. Figures [3] and 2] 
illustrate that incrtial-framc components of a tensor ex- 
pressed as functions of rotating-frame coordinates have 
some additional time-dependent mixing among the var- 
ious tensor-spherical-harmonic components. Therefore, 
filtering the inertial-frame components using the same 
algorithm used for single-coordinate frame evolutions is 
the wrong thing to do: it does not preserve all of the in- 
formation needed to determine a number of the highest- 
index coefficients in this case. This straightforward (and 
incorrect) implementation of spherical-harmonic filtering 
is the method used for the test shown in Fig. [2J with the 
result being unstable (in the sense that the constraints 
grow without bound) and probably non-convergent. The 
cure for this problem is also clear: transform spatial 
tensors to a rotating-frame tensor-sphcrical-harmonic ba- 
sis before filtering, then transform back to the inertial- 
frame basis afterwards. (Spacetime tensors, such as the 
four-metric ipatn are nrs t split into their spatial-tensor 
parts, e.g., V'tt; Vfe an d ^iji an d these parts are then 
filtered in this way.) Results using this new filtering al- 
gorithm are shown in Fig. [5] for the same initial data 
and dual-coordinate frames used in Fig. [21 These evo- 
lutions now appear to be stable and convergent on the 
needed timescale, except for mild (sublinear) power-law 
growth seen here only in the highest resolution cases. We 
suspect that this power-law growth may be due to ac- 
cumulated roundoff error, but we have not investigated 
this in detail because the growth is insignificant on the 
timescale needed for multiple-orbit binary black hole evo- 
lutions such as those described in Sec. [V] The new fil- 
tering method described here is needed to perform sta- 
ble evolutions of any rotating spacetime (even "linear" 
problems like evolutions of Minkowski spacetime) using 
the dual frame method. Additional filtering may also be 
needed in some circumstances to control other problems 
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FIG. 5: Constraint violations in dual-coordinate-frame evolu- 
tions of Schwarzschild, with comoving coordinates that rotate 
uniformly with angular velocity SI = 0.2 /M. These evolutions 
use the new rotating-frame tensor-spherical-harmonic filter- 
ing algorithm. The highest resolution case (dashed curve) 
uses more angular basis functions (i < 13 instead of I < 11) 
as a more stringent test of convergence. 

(spectral aliasing errors due to non-linearities for exam- 
ple), but no additional filtering was needed or used for 
any of the numerical evolutions presented in this paper. 



IV. BINARY BLACK HOLE TESTS 

The primary motivation for developing the dual- 
coordinate method describe in Sec. HT1 was to allow us to 
perform binary black-hole evolutions using coordinates 
that move with the holes. In this section we describe the 
first test of this new method for binary black-hole evolu- 
tions. We describe the binary black- hole initial data that 
we use, give a very brief description of the binary-black- 
hole specific features of our numerical methods, and then 
describe the results of this first test. 

The binary black-hole evolutions described here begin 
with initial data prepared with the methods described in 
Ref. [22J . These data represent two equal- mass corotating 
black holes in a quasi-stationary circular orbit. The data 
are obtained by solving the extended conformal thin- 
sandwich form of the initial value equations with the fol- 
lowing choices of freely specifiable data: the spatial con- 
formal metric is chosen to be flat, the trace of the extrin- 
sic curvature, its time derivative and the time derivative 
of the conformal metric are set to zero. In addition to de- 
termining the spatial metric and extrinsic curvature, this 
form of the initial value equations also determines a lapse 
and shift that produce relatively time-independent evo- 
lutions. We use the Neumann form of the lapse boundary 
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condition described in Ref. (22j. The binary black holes 
used for the tests described here have equal irreducible 
masses, M; 1T = 1.061536 (in our code coordinate units), 
total ADM energy M AD m = 2.100609 = 1.978834 M irr , 
and total ADM angular momentum Jadm = 4.3485850 = 
0.985502 M| DM . The coordinate separation of the cen- 
ters of these holes is set to 20 initially, and the initial spins 
of the holes are set to corotating values (i.e., no spin in 
the corotating frame). The initial orbital angular velocity 
of this binary is Q = 0.01418276 = 0.02979244/M ADM . 
This initial data set is from a family of publicly available 
initial data (23|. 

We represent the binary black holes considered here on 
a computational domain divided into 44 subdomains: 14 
spherical shells, 24 cylindrical shells, and 6 rectangular 
blocks. The various dynamical fields are expanded in the 
appropriate spectral basis functions for each subdomain. 
The evolutions described here are performed at two nu- 
merical resolutions using a total (over all subdomains) of 
260, 756 « 64 3 and 431, 566 « 76 3 collocation points re- 
spectively. Constraint preserving and physical boundary 
conditions are imposed at the outer boundary of our com- 
putational domain (initially at r = 280) as described in 
Ref. [f|, and the appropriate characteristic fields are ex- 
changed between subdomains at internal boundaries. A 
more detailed description of the numerical methods used 
for these binary black-hole evolutions will be included in 
a subsequent paper. Here we focus our attention on the 
dual-coordinate aspects of these evolutions. 

As our first binary-black-hole test of the dual- 
coordinate evolution method, we evolve the equal-mass 
circular-orbit binary described above using a rotating co- 
ordinate frame as described in Sec. IIIIl We set the an- 
gular velocity of this frame to = 0.02979244/M A dm, 
the initial orbital angular velocity of the binary. We 
track the position of the center of each black hole, 
[x c (t),y c (t), z c (t)], during this evolution by solving for 
the apparent horizon of each hole. We express the ap- 
parent horizon as a spherical- harmonic expansion [241 ] : 
-Rah = c + fJ2e, m ^rn(t)Yi m (6,ip), where c is the point 
about which we expand, f is a radial unit vector field cen- 
tered on c, and 9 and tp label the points on the apparent- 
horizon surface. The coefficients Ri m are obtained by 
minimizing the expansion at the Gauss-Legendre collo- 
cation points with a DFP-minimization [2g, Hf|. The 
center of the hole is definded as that c for which the 
£ = 1 components of the spherical harmonic expansion 
vanish: R\ m — 0. Initially the black holes are located at 
[x c (0),y c (0),z c (0)} = [±10,0,0]. The reflection symme- 
try of the initial data ensures that z c (t) = throughout 
the evolution, but x c (t) and y c (t) are free to drift. 

Figure [S] illustrates the motion of the center of one of 
the black holes, x c (t) and y c (t), with respect to the co- 
moving coordinates. We see that the center of the hole 
remains relatively fixed during the very first part of this 
evolution, but fairly quickly it begins to drift away from 
its initial location with x c (t) drifting more quickly than 
y c (t)- It appears that the two black holes get closer to- 
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FIG. 6: Evolution of the apparent horizon center, {x c ,y c ), 
of one of the black holes relative to a uniformly rotating co- 
ordinate frame. The minimum coordinate distance between 
the apparent horizon and the excision boundary becomes less 
than 0.003 at t — 45.6 Madm, and the evolution is terminated. 



gether due to the effects of gravitational radiation on the 
orbit. Our code automatically terminates an evolution 
when an apparent horizon gets too close (within 0.003 for 
this run) to an inner (excision) boundary of the computa- 
tional domain. If the apparent horizon were to cross this 
boundary then additional boundary conditions would be 
needed there and we do not know how to specify such 
boundary conditions in a physically meaningful way. 

Figure [7] shows the equatorial (x-y plane) cross sec- 
tions of the initial apparent horizon (solid curve), the 
final apparent horizon (dashed curve), and the location 
of the excision boundary (dotted curve) for one of the 
black holes from the evolution shown in Fig. [6] The 
apparent horizon has shifted to the left, confirming our 
diagnosis that the two black holes have gotten closer to- 
gether. Our code does not crash at the end of this test, 
and the evolution would have run (a little) longer had a 
smaller excision boundary been used. A smaller excision 
boundary results in higher truncation errors for a fixed 
numerical resolution, so making it smaller requires signif- 
icantly more computational resources. But even if cost 
were not an issue, making the excision boundary smaller 
would only delay the apparent horizon crossing the exci- 
sion boundary by a very brief time. A better solution is 
needed for this problem. 



V. HORIZON-TRACKING COORDINATES 

The binary black-hole test described in Sec. IIVI shows 
that a uniformly-rotating comoving coordinate frame is 
not adequate to keep the black holes centered on the ex- 
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FIG. 7: Moving coordinate frame representation of an equa- 
torial slice through the apparent horizon of one of the black 
holes for the evolution of Fig. [6] The evolution is terminated 
when the apparent horizon approaches the excision boundary 
(dotted curve) near the right side of this figure. 

cision boundaries of the computational domain. The mo- 
tions of the centers of the holes are too time-dependent 
for this simple approach to work effectively for very long. 
If we wish to construct coordinates that track the appar- 
ent horizons for many orbits, then a more flexible coordi- 
nate map is needed. Since the positions of the black holes 
are difficult (impossible) to predict a priori, some kind of 
feedback control system will also be needed to measure 
these positions and make the appropriate adjustments to 
the coordinate maps as the evolution proceeds. 

The motion of the black holes in the test of Sec. IfVI 
had a larger x c {t) component, so we begin by construct- 
ing a coordinate map that more accurately tracks this 
component of the black hole's position. Let us define a 
control parameter, Q x (t) = [x c (t) — x c (0)]/x c (0), that 
gives a dimensionless measures the ^-component of the 
position of the hole. The idea is to measure Q x (t) as the 
binary evolves, and then use this information to adjust 
the function a(t) in the time-dependent coordinate map, 



t 
x 

y 

z 



t. 



a(t ) [a; cos(ilt ) + y sin(f2t )] , 
a(t) \—x sin(f2t ) + y cos(f2t )] , 
a(t)z, 



(11) 
(12) 

(13) 
(14) 



in such a way that Q x (t) remains sufficiently small. This 
map applies the uniform rotation used in Sec. IIVI com- 
bined with a re-scaling of the spatial coordinates that can 
be used to keep the x-coordinate separation of the holes 
fixed. 

We borrow ideas from the literature on mathematical 
control theory [2?J to design a feedback control system 



for the map parameter a. The basic idea is to change a 
in such a way that Q x is driven back to its "equilibrium" 
value, Q x = 0, whenever it drifts away. The first step is 
to determine the response of the control parameter SQ X 
to a small change in the map parameter 5a. In this case 
the relationship is rather simple: SQ X — —Sa/a. For 
maps that are close to the identity, a = 1, this relation- 
ship is well approximated by 5Q X « —6a. We also want 
to make sure that any adjustment we make to the coordi- 
nate map is sufficiently smooth that it does not interfere 
with our ability to solve the Einstein equations. In par- 
ticular we need the coordinate map to be at least C 2 
so that the transformed dynamical fields u a are at least 
continuous [see Eqs. ©-©]■ Therefore our control sys- 
tem will be allowed to adjust only d 3 a/dt 3 freely, and a 
will be obtained by time integration to ensure sufficient 
smoothness. Following the usual procedure in control 
theory, we choose d 3 a/dt 3 in a way that is determined 
by the measured values of Q x and its derivatives: 



fa 
dt 3 



aQ x +f3 



dQ" 
dt 



+ 7 



d 2 Q x 



(15) 



where a, j3 and 7 are parameters that we are free to 
pick. We choose these parameters in such a way that the 
"closed-loop" equation, 



da 
It 3 



d 3 Q" 
dt 3 



= aQ x + [3 



dQ x 
dt 



•7- 



d 2 Q x 
dt 2 > 



(16) 



has solutions that all decay exponentially toward the 
desired equilibrium value Q x — > 0. We use the values 
a = A 3 , [3 = 3 A 2 and 7 = 3 A, that result in the following 
closed-loop equation: 



d3QX =-A 3 Q*-3X> dQX 



-3A 



d 2 Q 3 



(17) 



dt 3 *° dt dt 2 ' 

The most general solution to this equation is 

Q x (t) = {At 2 +Bt + C)e' xt , (18) 

for arbitrary constants A, B, and C. All of these solu- 
tions decay exponentially toward the desired equilibrium 
solution at a rate determined by the parameter A > 0. 

We use these ideas now to construct a feedback control 
system for a. We pick a set of control times ti at which 
the expression for a will be adjusted. The control interval 
At = ti + i — ti is chosen to be shorter than the timescale 
on which Q x drifts away from its equilibrium value. We 
choose the map parameter a(t) in the time interval ti < 
t < t i+ i to be: 



/ s / \ dai (t — ti) 2 d' 

dt 



ft, 



{t-Uf 



dt 2 



2 dt 2 



A 



dt 



(19) 



where the constants a,, dai/dt and d 2 ai/dt 2 are the val- 
ues taken from the map in the previous interval tj_i < 
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FIG. 8: Evolution of the apparent horizon center of one of 
the holes for a binary black-hole evolution using a feedback 
control system that adjusts the moving coordinate frame to 
control the parameter [x c — a; c (0)]/a; c (0). Compare to Fig. [6] 
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FIG. 9: Moving coordinate frame representation of an equa- 
torial slice through the apparent horizon of one of the black 
holes for the evolution of Fig. [8] Evolution is terminated when 
the apparent horizon approaches the excision boundary near 
the bottom of this figure. 



t < U evaluated at t = U, and the constants Qf, dQf /dt 
and d 2 Qf/dt 2 are the values measured at t = ti. This 
choice of a(t) guarantees that our comoving coordinate 
map is C 2 in time, C°° in space, and it enforces the 
closed- loop equation exactly at the control times t = U. 
We begin our evolutions at to — by setting the ini- 
tial conditions do = 1 and dao/dt = d 2 a /dt 2 = Qq = 
dQv/dt = d 2 Qfi/dt 2 = 0. To ensure that the discontinu- 
ities in d 3 a/dt 3 do not affect the convergence of our code, 
we require that the ti occur at fixed times (independent 
of the numerical timestep), and we choose the timesteps 
so that the U always coincide with the beginning of a full 
numerical timestep. 

We test this feedback control system by evolving the 
same binary black-hole initial data discussed in Sec. IIVI 
We set the value of the control damping parameter 
A w 0.3/Madm, the control interval At w IMadMi 
and the timescale used to evaluate average values of 
the derivatives of Q x in Eq. (02} to ~ 0.5Madm- Fig- 
ure [5] shows the motion of the center of one black hole 
during this test. We see that the control system effec- 
tively keeps the x-coordinate of the center of the hole 
Q x = [x c — x c (0)]/x c (0) within acceptable bounds. But 
this test ends at about t = 73.4 Madm when the other 
(uncontrolled) component of the black hole's position 
y c /x c grows too large. The new control system signif- 
icantly extends this binary evolution, but the binary has 
still only completed about 0.34 orbits. Figure [9] shows 
that the apparent horizon has drifted upward at the time 
this test ends. As the binary inspirals towards merger the 
angular velocity of the orbit increases, so the fixed an- 
gular velocity of the co-moving coordinates is no longer 
able to track the positions of the holes with sufficient 



accuracy. 

Fortunately, the dual-coordinate evolution method is 
extremely flexible and it is easy to construct a comov- 
ing coordinate map capable of tracking both the orbital 
rotation and the radial motion of the black holes. For 
example, the map 



t 
x 

II 

z 



t, (20) 

a(i ) [x cos <p{i) + y sin <p(i)] , (21) 

a(i)\—xsva.<p(i)+ycDS<p(ij\, (22) 

a(t)z (23) 



includes a time-dependent rotation angle ip(t) that can 
be used to adjust the angular velocities of the holes, in 
addition to the time-dependent conformal factor a(t). An 
additional feedback control system is then needed to ad- 
just the map parameter (p through measurements of the 
control parameter Q v = y c /x c . The basic control equa- 
tion for this parameter is SQ y = — 5 (p. So we construct a 
feedback control system for ip(t) that is completely anal- 
ogous to Eq. (fT§|) . In the time interval < t < ti we 
set 



<p(t) 



tpi + (t- 

(t-u) 3 



l) dt 



d 2 Qf 
dt 2 



{t-U) 2 d 2 tp t 
2 dt 2 
,dQ y 



+ X 



dt 



+ X 6 



,(24) 



where the constants ipi, dfi/dt and d 2 ipi/dt 2 are the val- 
ues taken from the map in the previous interval < 
t < ti evaluated at t = ti, and the constants Qf , dQf /dt 
and d 2 Qfjdt 2 are the values measured at t = ti. This 



10 



4x10 



2x10 



-2x10 



-4x10 




200 400 
t/M 



ADM 

FIG. 10: Evolution of the apparent horizon center of one black 
hole for a binary black-hole evolution using a feedback system 
that adjusts the moving coordinate frame to control both the 
parameters y c /x c and [x c — x c (0)]/x c (0). Compare to Figs. [6] 
and U 
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FIG. 11: Motion of the inertial frame coordinates of the center 
of the apparent horizon (x c ,y c ) (solid curve) of one of the 
black holes from the evolution of Fig. 1101 The dotted curve 
shows the equatorial cross section of the apparent horizon at 
t = 661.4 Madm- The binary system has completed about 4.6 
orbits. 



choice of (p(t) guarantees that our comoving coordinate 
map is C 2 , and it enforces the closed-loop equation peri- 
odically at the control times t — U. We begin our evolu- 
tions at tg = by setting the initial conditions dipo / dt = 
n and ip Q = d 2 ip /dt 2 = Q y = dQ y /dt = d 2 Q y /dt 2 = 0. 

We test this enhanced feedback control system, which 
controls both Q x and Q y , by evolving the same bi- 
nary black-hole initial data used for the evolutions in 
Figs. [H] and [5] We use the same control system param- 
eters for the Q x subsystem as those discussed above. 
For the Q y subsystem we use the value of the control 
damping parameter A « 0.5/A/adM; the control inter- 
val At « 0.5 Madm, and the timescale used to evaluate 
the average values of the derivatives of Q y in Eq. |24)) 
to ~ 0.2 Madm- Figure [TO] shows the control parameters 
Q x and Q y for this evolution. We see that both con- 
trol systems work extremely well, and the binary now 
evolves until t = 661.4 Madm before the evolution stops. 
Figure [TT] shows the evolution of the inertial- frame co- 
ordinates [x c (t ), y c (t )] of the center of one of the black 
holes, from which we see that this evolution has com- 
pleted about 4.6 orbits. Figure [12] shows the normalized 
constraint violations for this evolution, at two different 
numerical resolutions. The upper curve uses a total of 
260, 756 ~ 64 3 collocation points and the lower curve 
uses 431,566 « 76 3 collocation points. The constraints 
are normalized here by dividing by the initial value of the 
norm of the derivatives of the dynamical fields. Figure [P2l 
suggests that our numerical evolutions are convergent un- 
til just before the code terminates. 

Finally, Fig. [T3] shows the shape of the apparent hori- 
zon of one black hole (expressed in comoving coordinates) 



at t = (solid curve) and at the end of the evolution 
t = 661.4 Madm (dashed curve). The conformal factor 
a(t) used to keep the black holes at a fixed coordinate 
separation also causes the apparent horizons to expand 
in comoving coordinates. As the binary evolves, the ex- 
cision boundary (which has a fixed coordinate radius in 
the moving frame) moves deeper and deeper into the in- 
terior of the black hole. At fixed numerical resolution 
this leads to rapidly increasing truncation errors that in 
turn generate constraint violations. The dominant source 
of constraint violations seen at the end of the evolutions 
in Fig. [T2] comes from the region within the apparent 
horizons of the two black holes. This problem was par- 
tially corrected by moving the location of the excision 
boundaries at the time t = 554.1 Madm: for each hole, a 
portion of the unphysical interior grid was removed (by 
interpolating the dynamical fields onto a smaller grid), 
thus moving the excision boundary to the location of the 
outer dotted curve seen in Fig. [TBI This evolution stops 
in part because truncation errors grow rapidly as the ex- 
cision boundaries move deeper inside the holes. Another 
reason this evolution stops can also be seen in Fig. [13] 
the shape of each apparent horizon has become very dis- 
torted by tidal interaction with the companion hole. This 
distortion requires a significant increase in angular res- 
olution to represent it accurately; so at fixed numerical 
resolution, truncation errors quickly grow and these cause 
rapidly growing constraint violations seen near the end 
of the evolutions in Fig. [TJ] 
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VI. DISCUSSION 



— t=0 




I I I I I I I I I I I I I 

8 9 10 11 12 

x 

FIG. 13: Equatorial cross sections, in comoving coordinates, 
of the apparent horizon of one hole for the evolution of 
Fig. 1101 The horizon is shown both at the beginning and 
at the end of the evolution. The smaller dotted circle (just 
inside the initial horizon) shows the excision boundary used 
for t < 554.1 Madm- The outer dotted circle shows the exci- 
sion boundary used for t > 554.1 Madm- 



We have shown in Sec. IIIII that the dual-coordinate 
frame evolution method introduced in Sec. [TT] is capa- 
ble of solving the rotating-frame instability problem. A 
more sophisticated application of this method has also al- 
lowed us to construct comoving coordinates in Sec.lVlthat 
track the apparent horizons of a binary black hole sys- 
tem through a feedback control system. These horizon- 
tracking coordinates allowed us to evolve a binary black- 
hole system stably and accurately for about 4.6 orbits. 
Our expectation is that the great flexibility of the dual- 
coordinate method will make it useful in many other ap- 
plications. We imagine, for example, that coordinate 
maps could be constructed that keep the shape of the 
apparent horizons spherical, and their locations close to 
the excision boundary. Such maps should solve the prob- 
lems that prevent the binary black-hole evolution shown 
here from proceeding all the way to merger. We also ex- 
pect that many applications are possible that we are not 
able to anticipate at this point. 



Appendix 

In this appendix, we describe some of our efforts to 
understand and cure the instability that occurs when 
solving the generalized harmonic system using a single 
rotating coordinate frame. We found the instability, il- 
lustrated in Fig. [TJ to occur even for evolutions of fiat 
spacetime in a rotating frame. Since Minkowski space 
is simpler than Schwarzschild, much of our analysis was 



focused on understanding the flat-space case. 

At the analytical level, it is easy to see that in 
Minkowski space stability in non-rotating inertial coordi- 
nates is equivalent to stability in a single rotating coor- 
dinate frame. The basic argument is as follows. Con- 
sider two coordinate systems x a and x a related by a 
uniform rotation, Eqs. (|7|)- (|10p . The evolution equa- 
tions for linear perturbations of the metric 5ip a b about 
a flat Minkowski background may be written in an ab- 
stract way as E a b[Sip] — 0. We assume that both the 
fields and these evolution equations transform as tensors, 
i.e., 5ip ab = A a a Aif5ip ah and E ab = A a a A b l E al , where 
A a a = dx a /dx a . This implies that if 5ip db - is a solution in 
the inertial frame then 5ip ab = A a a A b b 5ip a i is a solution 
in the rotating frame. Suppose that 5ip ab ~ e st e m ^ is an 
eigenmode in the inertial frame. Then 8ip a b ~ e s t e lmv 
is a solution in the rotating frame with s' = s — im£l, 
so that Re(s') = Re(s). Thus, linear stability in the two 
frames is equivalent. We have done a complete linear 
stability analysis in the inertial frame, and the solutions 
are (of course) stable. Nevertheless, our generalized har- 
monic evolution code is stable for nearly flat evolutions 
in non-rotating coordinates, but is extremely unstable for 
solutions in a rotating frame. 

To understand this apparent contradiction, let us now 
consider the second-order form of the vacuum generalized 
harmonic evolution equations in more detail, 

E ab [i)] = RabbP] - V (a C 6) [<0] = 0, (25) 

where ip ab is the four-metric, V a its covariant derivative, 
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R ab the Ricci tensor, and C a is the harmonic constraint 
defined in Eq. ©. Equation (|25[) transforms as a tensor 
if and only if C a transforms as a vector. For this to be 
true, we have to demand a rather special transformation 
law for the gauge source function H a : from the known 
transformation of the Christoffcl symbols and Eq. ([6]) we 
deduce 

H a = H a - M hc t d bc , (26) 

where H a = A a a H- a and f Q bc = (A^ 1 ) a ad c A b a . Naively, 
in order to evolve Schwarzschild spacetime in a rotating 
frame, say, one would just set H a to be the — T a of the 
exact solution evaluated in that frame and freeze it in 
time. In contrast, Eq. ((26)) tells us that H a should depend 
on the evolved metric ip ab . We find that freezing the 
quantity H a in the rotating frame and then obtaining 
H a from Eq. (|26[) using the evolved metric ?p ab makes our 
rotating-frame Minkowski and Schwarzschild evolutions 
much more (but not completely) stable. 

Our numerical code solves a first-order representation 
of the generalized harmonic evolution equations. To un- 
derstand the remaining problems we must look at the 
first-order reduction in more detail. The first-order for- 
mulation introduces the new dynamical fields 

U ab = -t c d c ip ab , $ fca6 = d k ip ab . (27) 

Unlike the metric ip ab , these new fields do not transform 
as tensors, as can be appreciated from Eqs. (|3|)-(j4|). Thus 
the simple analytical argument given above for the stabil- 
ity of rotating flat space does not apply. We can modify 
the standard first-order form of the equations by chang- 
ing the definitions of these new fields slightly. We replace 
Eq. ((271) witn 

IT Qb = -t C V c tlj ab , $fcab = Vfc^ab, (28) 

where V a denotes the covariant derivative associated 
with the connection T a bc introduced in Eq. (f26|) . By def- 
inition, V a reduces to d a in the inertial frame, but not 
in the rotating frame. The modified variables in Eq. (|28|) 
by definition transform as tensors under rotations of the 
coordinates. To make this new first-order formulation 
truly covariant we must also ensure that the evolution 
equations (and boundary conditions) transform as ten- 
sors. To do this, we start with the equations written in 
the inertial frame (see e.g., Eqs. [35] -[37] of Ref . 0) and 
formally replace all partial derivatives d a with the co- 
variant derivatives V a . The rotating frame solutions of 



these modified evolution equations should then have the 
same stability properties as the solutions obtained in the 
inertial frame. We have used this system to evolve both 
Minkowski and Schwarzschild spacetimes, and found the 
numerical solutions to be greatly improved over the un- 
modified system. But the results are still not quite stable. 

The next problem that we discovered turned out to 
be numerical: we needed to modify the standard tensor- 
spherical-harmonic filtering algorithm described for ex- 
ample in Sec. IIIII The transformation that relates the 
rotating frame components of tensors with non-rotating 
frame components mixes together the time components 
with some spatial components (e.g., tptt depends on tpn, 
i\)t<p and t/'w)- This causes scalar, vector, and second- 
rank tensor spherical harmonics to be mixed together in 
a non- linear way. To disentangle these parts, it is nec- 
essary first to transform the time components of these 
fields to the inertial frame, filter them there, and then 
transform them back to the rotating frame. In this way, 
we remove the mixing of the different tensor spherical 
harmonics that occurs in a rotating frame. Using this 
filtering method (which turns out to be equivalent to 
the filtering used for our dual-coordinate evolutions de- 
scribed in Sec. IHIj) produces Minkowski-space evolutions 
that appear to be completely stable, while evolutions of 
Schwarzschild though improved are still not quite stable. 

The combination of the modifications described above 
significantly improves the stability of our evolutions of 
flat space and single black-hole spacetimes in rotating 
coordinates. However, these methods did not produce 
robustly stable evolutions over the entire range of do- 
main sizes and rotating frame angular velocities needed 
for binary black hole simulations. We stopped pursuing 
this approach after the dual-frame method was found to 
be so stable, and so flexible that it allowed us to solve 
our moving frame excision problem as well. 
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